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Abstract Using a simple color selection based on B-, z- and if-band photometry, BzK — 
{z — K)ab — {B — z)ab > —0.2, we picked out 52 star-forming galaxies at 1.4 ^ z 2.5 
(sBzKs) from a K-bwd selected sample (AVcga < 22.0) in an area of ^ 5.5 arcmin^ of 
the Hubble Ultra Deep Field (UDF). We develop a new photometric redshift method, and 
the error in our photometric redshifts is less than 0.02(1 + z). From the photometric redshift 
distribution, we find the BzK color criterion can be used to select star-forming galaxies at 
1.4 < z <, 2.5 with Kycga < 22.0. Down to Xvcga < 22.0, the number counts of sBzKs 
increase linearly with the K magnitude; the sBzKs are strongly clustered, and most of them 
have irregular morphologies on the ACS images. They have a median reddening of E{B — 
V) ^ 0.28, an average star formation rate of ~ 36 Mq yr^^ and a typical stellar mass of 
~ 1O^°M0. The UV criterion for the galaxies at z ^ 2 can select most of the faint sBzKs in 
the UDF, but it does not work well for bright, massive, highly-reddened, actively star-forming 
galaxies. 

Key words: galaxies: evolution — galaxies: high-redshift — galaxies: photometry — cos- 
mology: observations 



1 INTRODUCTION 

A number of observations suggest that the era of z 2 is important in galaxy evolution for various reasons: 
the cosmic star formation rate begins to drop at z 1 — 2 from a flat plateau at higher redshifts; the 
morphological type mix of field galaxies changes remarkably at z ~ 1 — 2; the number density of QSOs has 
a peak at z - 2 (Dickinson et al. 2003; Fontana et al. 2003; Steidel et al. 2004; Kong et al. 2006; Richards 
et al. 2006). However, galaxies in the redshift interval 1.4 ^ z ^ 2.5 have been called the "redshift desert" 
by some researchers because, to date, we lack the tools to pick out dim galaxies in that range. 

Recently, Steidel's group has extended the UV technique for selecting LBGs (Lyman Break Galaxies) 
to z < 3 using a UnGRs color-color diagram which isolates star-forming galaxies ai z ^ 2 (Erb et al. 2003; 
Steidel et al. 2004; Adelberger et al. 2004). However, star-forming galaxies can be selected as LBGs only if 
they are UV bright (i.e. actively star forming) and not heavily reddened by dust. Using the highly complete 
spectroscopic redshift database of the K20 survey, Daddi et al. (2004) introduced a new criterion, BzK = 
(z — K)p^-Q — {B — z)ab > —0.2, in the B-, z- and if -band photometry, for obtaining a virtually complete 
sample of star-forming galaxies in the redshift range 1.4 ^ z ^ 2.5 (for convenience, we use the term sBzKs 
for the galaxies selected in this way). This criterion is reddening independent for star-forming galaxies in 
the selected redshift range, thus can be used to select the reddest dust-extinguished galaxies. This should 
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allow a relatively unbiased selection of star-forming galaxies at z ~ 2 within the magnitude hmit of the 
sample studied. 

Based on BRIzJK photometry obtained by combining Subaru optical and ESO near-IR data over two 
separate fields and using the BzK selection technique, Kong et al. (2006) has obtained complete samples 
of ~ 500 candidate sBzKs which were identified over an area of ~ 920 arcmin^ to -fCvoga = 19.2, of which 
320 arcmin^ are complete to Kvoga = 20. Using HST/ACS-/ (F814W), Subaru-B^g'^Vz', CFHT-wi and 
KPNO/CTIO-Xs data, Renzini et al. (2008) has selected ~ 1300 sBzKs over an area of - 6600 arcmin^ to 
K^ega. = 19-2 in the COSMOS survey (Scoville et al. 2007). The surface density, number counts, clustering, 
reddening, star formation rates (SFRs) and stellar masses (M*) of these bright sBzKs (ii'vega < 20 or 
Kyega < 19.2) Were analyzed in these papers. To select/am? sBzKs and study their properties, deep infrared 
images are required. The Hubble Ultra Deep Field (UDF) provides us with the deepest view to date of the 
visible universe. In this paper, we will use the space-based and ground-based UDF data to select a faint 
sBzKs sample (down to iv'vega = 22) to study the properties of faint sBzKs, and compare the differences 
between faint and bright sBzKs. 

This paper is organized as follows. Section 2 describes the UDF observations and the method for ob- 
taining the photometric catalog. Section 3 describes the calculation and calibration of photometric redshifts. 
Section 4 presents the selection, morphologies, number counts, clustering and physical properties of sBzKs. 
Section 5 compares the samples selected with the BzK technique to those UV-selected galaxies at z ~ 2. 
Finally, a brief summary is presented in Section 6. For the sake of comparison with previous works, magni- 
tudes and colors in both AB and Vega systems had to be used. 

2 UDF OVERVIEW 

The UDF (RA=03^32'"39^0, Dec=-27°47'29'.'l, J2000) is located within one of the best studied areas of 
the sky: the Chandra Deep Field South (CDF-S or GOODS-South region). With a total of 544 orbits, it is 
one of the largest time allocations with HST, and indeed the filter coverage, depth, and exquisite quality of 
the UDF ACS and NICMOS images provide an unprecedented data set for the study of galaxy evolution. 
The field has been imaged by a large number of telescopes at a variety of wavelengths. In this paper, three 
sources of imaging data and the VLT/FORS2 spectroscopic data for this field are used. 

2.1 Observations 

The UDF is a public HST survey of a single Advanced Camera Survey (ACS) wide field camera (WFC) 
field (11.5 arcmin^) in 4 broad-band filters: F435W (B), F606W (V), F775 W (i) and F850LP (z) (Beckwith 
et al. 2006). For our analysis we use the reduced UDF data vl .0 made public by the Space Telescope Science 
Institute on 09 March 2004. The lO-cr limiting magnitudes at i?-band, V-band, i-band and z-band are 28.7, 
29.0, 29.0 and 28.4, respectively, in an aperture of 0.2 arcsec^. 

We also use the HST/NICMOS (Near-Infrared Camera and Multi-Object Spectrometer) UDF photom- 
etry data for our analysis, which were taken with NICMOS under the HST Cycle 12 Treasury Program in 
two bands, Fl lOW and F160W, which roughly correspond to J and H. Due to the small field of NICMOS 
camera, the NICMOS UDF only covers a part (5.76 arcmin^) of the optical UDF. The data reduction and 
photometry were introduced by Thompson et al. (2005). The 5-a limiting AB magnitude is 27.7 at 1.1 and 
1.6 /im in a 0.6" diameter aperture. 

Ground-based near-infrared images (JHK) of the UDF were taken as part of the GOODS survey with 
ISAAC on the VUT. These data were obtained as part of the ESO Large Programme LP168.A-0485(A), 
and GOODS/ISAAC Data Release Version 1.5 was used for our analysis. This data release includes 24 
fully reduced VLT/ISAAC fields in J and K bands covering 159.1 and 159.7 arcmin^ of the GOODS 
region, respectively, and 19 fields in H band covering 126.7 arcmin^. The ISAAC data were reduced using 
an improved version of the ESO/MVM image processing pipeline. The 5-a limiting AB magnitudes at 
J-band, i7-band and if -band are 25.3, 24.8 and 24.4, in a 2.0" diameter aperture. 

A catalog of objects whose redshifts are known from spectroscopic observations is very important to 
the construction of a photometric redshift code. For this reason, the GOODS multi-object, optical to NIR 
spectroscopy in CDF-S was used for our analysis (Version 2.0, 24 December 2005), with the FORS2 on the 
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VLT. We obtained 1204 spectra of 930 individual targets, providing in total 943 redshift measurements with 
quality flag A, B or C (A=solid redshift, B=likely redshift, C=potential redshift). Of the targets 725 have 
an assigned redshift with quality flag A, B or C. A full description of the survey can be found in Vanzella 
et al. (2005). Forty-eight of the targets are located in the UDF, of which 23 have quality flag "A" and 10 
have quality flag "B". 

2.2 Image Resampling 

All of these data have been released in fully processed form, and no additional processing is necessary. 
However, the images in different data sets have different scales and sizes, so they must be resampled in 
order to put on the same astrometric grid. The resampling was done with IRAF's geomap and geotran 
tasks. The NICMOS H image was used as the base image to which all the other images were matched. 
The NICMOS J image need no remapping as it was already on the same scale. The optical images from 
the HST/ACS and the near-infrared images from VLT/ISAAC were remapped. This changed their scales 
from 0.03 and 0.15 to 0.09 " pixel^^. The resampling does not introduce any appreciable shifts in either the 
position or flux. 

Figure [T] shows a composite pseudo-color image of the UDF. From this figure, we find that the full 
drizzled image does not have a uniform integration time over the image. In particular, the edges of the 
HST/NICMOS image have only one integration, as compared to the average 16 integrations for the interior 
of the image. Therefore, the area, as discussed in this paper, was reduced from HST/NICMOS 's 5.76 to 
5.50 arcmin^, as is outlined near the center of the image. 

2.3 Photometry 

Source extraction in the science image was performed with the source extraction program SExtractor 
version 2.5 (Bertin & Arnout 1996) in the dual image and rms image mode. The source extraction includes 
four HST/ACS images, three VLT/ISAAC images, as well as two HST/NICMOS images. The NICMOS 
J+H image was used as the detection image for our catalog. The source extraction parameters were similar 
to those used in the NICMOS treasury version 2.0 catalog, except for 2" diameter aperture used for the 
aperture magnitudes. 620 sources were detected in the central 5.5 arcmin'^ area down to Hab — 26.0 mag. 
For most of the objects, our Kron automatic aperture magnitudes (MAG_AUTO) for BVizJH are close to 
those in Thompson et al. (2005), with differences generally less than 0.05 mag. 

The total magnitude was then defined as the brighter of the Kron automatic aperture magnitude and the 
corrected aperture magnitude. The colors were measured using the 2" diameter aperture magnitudes. All 
magnitudes were corrected for Galactic extinction. As — 0.033, taken from Schlegel et al. (1998), using 
the empirical selective extinction function of Cardelli et al. (1989). Compared to the optical selection, the 
near-lR selection (in particular in the K band) offers several advantages (Cimatti et al. 2002). Therefore, 
we select objects to Kycga < 22.0 over a total sky area of 5.5 arcmin^ in the UDF, and 210 objects were 
included in our final catalog. Simulations of point sources show that in all the area the completeness is 
above 90% at this iC-band level. 

Figure |2] shows a comparison of the /-C-band number counts in the UDF survey with a compilation 
of counts published in the literature. No corrections for incompleteness were applied to our data, and we 
excluded stars using the same method as that in Kong et al. (2006). The red-, black-, green- and blue-filled 
squares correspond to the counts of field galaxies in the UDF, COSMOS (Renzini et al. 2008), Daddi-F and 
Deep3a-F (Kong et al. 2006), respectively. As shown in the figure, our number counts in different fields are 
in good agreement with those of the previous surveys. 

3 PHOTOMETRIC REDSHIFTS 

Photometric redshifts (hereafter Zphot) were calculated for all the Xvcga ^ 22 objects in the UDF. We take 
the public code hyperz, which uses a usual template-fitting and a minimization method (Bolzonella 
et al. 2000). 

We use a stellar population synthesis model (Kodama & Arimoto 1997, hereafter KA97) to make tem- 
plate SEDs. KA97 was successfully used to obtain photometric redshifts of low redshift galaxies (Kodama, 
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Fig. 1 Composite pseudo-color image of the UDF. The RGB colors are assigned to H-,z-, and 
B-band images from HST. The outlined (white) region near the center of the image is the field 
where the Ji7-band images have a high signal-to-noise ratio (5.5 arcmin^). 



Bell & Bower 1999) and high redshift galaxies (Furusawa et al. 2000). The template SEDs consisted of the 
spectra of pure disks, pure bulges and composites obtained by interpolation. In total, our basic template set 
consisted of 303 SEDs (more description in a forthcoming paper). The spectral energy distribution derived 
from the observed magnitudes of each object was compared to each template spectrum (redshift from 0.0 
to 5.0 at steps of 0.05; Ay from 0.0 to 4.5 at steps of 0.05 and Calzetti et al. 's internal reddening law) in 
turn. The weighted mean redshift from hyperz was then calculated for each object. 

Figure [3] shows a comparison of the photometric redshifts with the 33 spectroscopic redshifts from 
VLT/FOS2 spectroscopy (Sect. [ZTT i. Filled squares show the objects with "solid redshifts" (quality flag A) 
while open squares show the objects with less secure spectroscopic redshift identifications (quality flag B, 
"likely redshift"). Taking into account only the objects with solid redshifts, the error in the photometric 
redshifts is ct^ — 0.003(1 + z) with no catastrophic failures, but taking into account all the 33 objects, the 
error would increase to az = 0.02(1 + z). 
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Fig. 2 if-band differential number counts for field 
galaxies, EROs and sBzKs (see Sect. 4) in the UDF, 
compared with a compilation of results taken from 
various sources. Filled squares, triangles and open 
squares represent the number counts of if-selected 
field galaxies, EROs and sBzKs, respectively. 



Fig. 3 Comparison between photometric (Zphot) and 
spectroscopic redshifts (zspec) for the UDF spectro- 
scopic sample. Filled squares and open squares repre- 
sent respectively objects with solid redshifts (quality 
flag A) and likely redshifts (quality flag B). 



4 sBzKs IN THE UDF 



In this section, we select sBzKs in the UDF using the multicolor catalog based on the NIR fT-band image. 
Then we will discuss the number counts, clustering, morphologies, reddening, SFR and stellar masses of 
faint sBzKs in the UDF. 



4.1 sBzKs Selection 

The UDF z- and /f -band filters are the same as those used in Daddi et al. (2004) while the i3-band filters 
are different in the UDF (F435W) and Daddi etal.(2004) (Bessel B). Using the Pickles (1998) stellar 
spectra, we find the correction term to the B-band is very small (< 0.05 in all cases). So we applied the 
BzK selection criterion from Daddi et al. (2004) without change to the sBzKs selection in the UDF without 
additional correction. 

Figure St shows the BzK color diagram of ii'-selected objects in the UDF. Using the color criterion 
from Daddi et al. (2004), BzK = {z ~ K)ab - {B - z)ab > -0.2, 52 galaxies with Xvcga < 22 are 
selected in the UDF as sBzKs, which occupy a narrow range on the left of the solid line in Figure 4a. To 
^Vega < 22, the surface density of sBzKs is ^ 9.45 ± 1.05 arcmin"^. We also selected EROs in the UDF 
with the criterion i ~ K > 2.4 (AB mags), corresponding closely to i? — A' > 5 (Vega mags). In the UDF 
to Kycga < 22 24 EROs are selected and plotted in Figure 4a with open circles. 

FigurellJ) shows the photometric redshift histogram of 52 -ft^vcga ^ 22 UDF galaxies selected with the 
BzK criterion. The average redshift and the peak of the redshift distribution of sBzKs are 2;phot ^ 1-8. Only 
about 10% of them (6 objects) have Zphot < 1-4, three of these six objects are EROs, and the other three are 
just marginally inside the BzK selection region. The photometric redshift distribution supports the view 
that the BzK criterion is valid for picking out star-forming galaxies at 1.4 <, z < 2.5 with ATvcga — 22.0. 
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Fig. 4 a) The {z — K)ab vs {B — z)ab color-color diagram for the galaxies in the UDF. 
The diagonal solid line defines the region BzK = {z — K)ab — {B — z)ab > —0.2 that 
is efficient for isolating z > 1.4 star forming galaxies (sBzKs). The dashed line separates the 
regions occupied by stars and galaxies. Solid circles show objects classified as stars, having 
(z — A')ab — 0.3 (-B — z)ab < —0.5; open stars show stellar objects from the K20 survey 
(Daddi et al. 2004); squares represent sBzKs; circles represent galaxies with {i — K)ab > 2.40 
(EROs); solid stars correspond to stellar objects given in Pickles (1998). b) Photometric redshift 
histogram of 52 Kyoga. ~ 22 UDF galaxies selected with the BzK criterion (sBzKs). 



4.2 Morphologies of sBzKs 



The morphology of a galaxy reflects its dynamical history and evolution. For example, different Hubble 
types are associated with different star formation histories and different patterns of motion of stars and gas. 
Morphological studies of distant objects have become possible now with the high-resolution capabilities 
of HST/ACS. The ACS imaging provides a fundamental complement of investigating the nature of sBzKs 
and elucidating their evolutionary status. Figure|5]shows the color images for 48 sBzKs (to show the image 
well, we randomly selected 48 sBzKs from all 52 sBzKs) in the UDF. The B-, i- and z-band images were 
colored red, green, and blue, respectively. 

From this figure, we find that although the sBzKs show a range of morphologies from compact protob- 
ulges to disturbed major-mergers, most of them show irregular/merging-like morphologies. The bottom row 
in Figure 5 shows objects which satisfy both the BzK > —0.2 criterion (sBzKs) and the i — K > 2.4 crite- 
rion (EROs). More detailed analysis, including concentration (C), asymmetry (A), dumpiness (5), the Gini 
coefficient (G) and M20 for the sBzKs and comparison with other high redshift galaxies will be presented 
in our forthcoming paper. 
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Fig. 5 Composite pseudo-color images (3" x 3") of 48/52 sBzKs in the UDF. The RGB colors 
are assigned to z-,i- and S-band images from HST/ACS. Most of sBzKs show irregular/merging- 
like morphologies. 



4.3 Number Counts of sBzKs 

The differential number counts of sBzKs and EROs in the UDF are plotted in Figure |2l Open squares and 
solid triangles with solid lines in Figure 2 show the number counts for the sBzKs and EROs in the UDF, 
respectively. 

From Figure 2 some characteristics are distinguished: 1) The fraction of sBzKs increases very steeply 
towards fainter magnitudes, and the counts of sBzKs have roughly the same slope at all K-band magnitudes. 
2) The slope of the number counts of EROs is variable, being steeper at bright magnitudes and flattening 
out towards faint magnitudes. A break in the counts is present at iiTvcga ~ 18.0, very similar to the break 
that has been shown in the previous works. 3) The counts of field galaxies, EROs and sBzKs in the UDF 
are almost identical to those in COSMOS, Daddi-F, and Deep3a-F to their limits of Ky^ga ~ 20. 

4.4 Clustering of sBzKs 

Measurement of the clustering of galaxies provides an additional tool for studying the evolution of galaxies 
and the formation of structures. In this section we estimate the angular correlation of the general galaxy 
population as well as that of sBzKs. In order to measure the angular correlation function of various galaxy 
samples, we apply the Landy-Szalay technique (Landy & Szalay 1993; Kerscher et al. 2000). A fixed slope 
5 — 0.8 is assumed for the two-point correlation function. 
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In Figure |6h, the bias-corrected two-point correlation functions w{d) of field galaxies are shown as 
squares. The dashed line shows the power-law correlation function given by a least squares fitting. We 
clearly have detect a positive correlation for the field galaxies with an angular dependence broadly consistent 
with the adopted slope S — 0.8. The derived clustering amplitudes (where A is the amplitude of the true 
angular correlation at 1 degree) are A = 1.54, 1.46, 1.40, 1.30 and 1.26 for Ky^^^ < 20.0, 20.5, 21.0, 21.5, 
22.0, respectively. 

We estimate the clustering properties of the UDF sBzKs (the triangles in Figure 6a). The dotted lines 
show the power-law correlation function given by a least squares fitting to the measured values. A strong 
clustering of sBzKs is indeed present at all scales, with amplitudes about 6 times higher than those of the 
field population at the same i^^vcga limits, in agreement with the previous findings (Kong et al. 2006). The 
derived amplitudes are A = 8.07, 7.06, 6.30, 5.59 and 4.68 for Ky^g^ < 20.0, 20.5, 21.0, 21.5, 22.0, 
respectively. 

Figure 6b summarizes the clustering measurements for the two populations examined (field galaxies 
and sBzKs), as a function of the X-band limiting magnitudes of the samples. Clear trends with the K 
magnitudes are present for all samples, that the fainter galaxies apparently have lower intrinsic (real space) 
clustering, consistent with the fact that objects with fainter K magnitudes are less massive, or have wider 
redshift distributions, or both. We compare our results on the UDF with those previously reported by Renzini 
et al. (2008) for the COSMOS and Kong et al. (2006) for the Daddi-F and Deep3a-F, and find they all show 
a smooth decline in amplitude with X-band magnitudes. However, the decline is not less steep in the range 
18 < A'vcga < 20. Furthermore, the clustering amplitudes in these fields are slightly different. The field-to- 
field variation may be one of the possible reasons for this discrepancy, depending on the survey geometry, 
surface density and clustering properties. 

4.5 Physical Quantities of sBzKs 

In this subsection, physical properties of sBzKs such as color excess E{B — V), SFRs and stellar masses, 
are derived on the basis of the present photometric data and the recipes calibrated in Daddi et al. (2004). 

Following Daddi et al. (2004), the reddening of sBzKs can be estimated by the B — z colors, providing 
a measure of the UV slope. The Daddi et al. 's recipe is consistent with the recipes by Meurer et al. (1999) 
and Kong et al. (2004) for the Calzetti et al. (2000) extinction law, based on the UV continuum slope. The 
E{B — V) histogram of sBzKs in the UDF is shown in Figure|7^ (solid line). The median reddening for the 
i^Voga < 22 sBzKs is estimated to be E{B -V) = 0.28, and 18% of sBzKs have E{B ~V)> 0.4, where 
the UV-based criterion of Steidel et al. (2004) would fail to select the galaxies at 2 ~ 2. For comparison, we 
show the E{B - V) histogram of sBzKs in the COSMOS (ifvcga < 19.2) in Figure|7}i (dotted line), too. 
The median reddening for the Xvcga < 19.2 sBzKs in the COSMOS is estimated to be E{B -V) = 0.42, 
and 56% of them have E{B — V)> 0.4. The difference is due to the different Hmiting K magnitudes in 
the two fields. 

Knowing the reddening, the reddening corrected B-band flux is used to estimate the 1500 A rest-frame 
luminosity, assuming an average redshift of 1.8, which can then be translated into SFR on the basis of the 
Bruzual & Chariot (2003) models. Daddi et al. (2005) showed that SFRs derived in this way were consistent 
with radio and far-IR based estimates for the average sBzKs. The SFR histograms of sBzKs in the UDF and 
COSMOS are shown in Figure 7b as solid line and dotted line, respectively. About 16% of the sBzKs in 
the UDF have SFR> 70 Mq yr"\ and the median SFR is about 36 Mq yr"^ However, about 98% of flie 
sBzKs in the COSMOS have SFR>70 Mq yr^^, and the median SFR is 430 Mq yr^^. 

Following Daddi et al. (2004), the stellar masses of the sBzKs are estimated by an empirical relation 
between their observed A' -band total magnitudes and z — K colors. The histograms for the stellar masses of 
sBzKs derived in this way are shown in Figure 7c. Almost all of the sBzKs in the UDF have Af* < IQ^^Mq 
and the median stellar mass is 1.4 x IO^^Mq; in the COSMOS, - 85% of the sBzKs have AU > IO^Mq 
and the median stellar mass is ^ 1.5x IO^^Mq. The higher masses for the sBzKs in the COSMOS compared 
to the UDF result from the shallower X-band limit. 
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Fig. 6 a)Observed, bias-corrected two-point correlations of field galaxies (squares) and sBzKs 
(triangles). The error bars on the direct estimated values are Icr errors. The lines show the power- 
law fitted to the w{6). b) Angular clustering amplitudes of field galaxies and sBzKs in different 
fields are shown as a function of A'-band limiting magnitudes. 



5 sBzKs- VS. UV-SELECTED GALAXIES AT Z - 2 

Recently, the UV technique for selecting LBGs has been extended to z < 3 using a UnGR color-color dia- 
gram, which isolates star-forming galaxies at 1.4 < z < 2.5 (Erb et al. 2003; Steidel et al. 2004; Adelberger 
et al. 2004), a redshift range matching that of the BzK selection. To compare these different color criteria 
(BzK-criterion and UV-criterion) for selecting star-forming galaxies at 1.4 < 2; < 2.5, we show the sBzKs 
in the UDF and COSMOS on the UnGR color-color diagram in Figured] 

As no UnGR photometry is available to us for the sBzKs in the UDF region, synthetic UnGR magni- 
tudes have been derived from the KA97 models which provided the best fits to the observed ACS-BViz, 
NIC- Ji? and ISAAC-JHK SEDs. Figure 8 shows the resulting synthetic (G - R) vs. (C/„ - G) colors. 
Among the 52 Xvcga < 22 sBzKs in the UDF, 39 (75%) would be selected by the Adelberger et al. (2004) 
criterion for the galaxies at 2: ~ 2, five would be selected by the Steidel et al. (2003) criterion for the galax- 
ies at z ~ 3, only eight i^voga < 22 sBzKs fail to be selected by the UV criterion (these are highly dust 
obscured galaxies). 

InFigure[8] we also plot the COSMOS sBzKs on the synthetic (G-i?) vs. (C/„-G) diagram. Synthetic 
UnGR magnitudes have been derived from the KA97 models that provided the best fits to the observed 
CFHT-ui, Suhaiu-BVg'r'i'z', SDSS-u' g' r' i' z' and KPNO/CTIO-if s SEDs. Among - 1300 K < 19.2 
sBzKs in the COSMOS, only 450 (35%) would be selected by the Adelberger et al. (2004) criterion while 
about 65% Kycga < 19.2 sBzKs fail to be selected by the UV criterion. 
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Fig. 7 Panels (a) - (c) show, respectively, the histograms of reddening, star formation rate and 
stellar mass. Solid histograms refer to sBzKs in the UDF (numeric label on the left), dotted 
histograms, to the COSMOS objects (numeric labels on the right). 




Fig. 8 UnGR color-color diagram for sBzKs in the UDF (squares) and COSMOS (stars). The 
UnGR colors were derived from SED fitting. The solid line shows color regions defined for the 
UV identification at z ~ 2 and z ~ 3. 



From Figure 8 we find that most of the faint sBzKs can be selected by the UV criterion for the galaxies at 
z ~ 2, but it does not work well for the brighter sBzKs objects. One possible reason is that the UV-selection 
requires the UV continuum to be relatively flat, thus limits the overall dust extinction to E{B — V) ^ 0.4 
(Adelberger & Steidel 2000), but many of the bright sBzKs at z > 1.4 are more reddened. As we discussed 
in Section 4.5, a large proportion of highly reddened galaxies are actively star-forming and have large 
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stellar masses. Therefore, a large fraction (as high as ~ 65%) of the K < 19.2 galaxies at z > 1.4 fail to 
be selected by the UV criterion of identifying galaxies at 2; ~ 2. As a consequence, a significant part of the 
SFR density and stellar mass density is missed. 

6 SUMMARY AND CONCLUSIONS 

In this paper, we investigate the properties of star-forming galaxies at lA ^ z ^ 2.5 (sBzKs) in the 
HST Ultra Deep Field. Using images from ACS-BViz, NIC- JH, ISAAC- JHK and spectroscopy from 
VUr/FOS2, we analyze the photometric redshifts, number counts, clustering, E{B — V), SFRs, stellar 
masses and morphologies of sBzKs in the UDF. The main results can be summarized as follows. 

1 . Based on the KA97 templates and a standard minimizing method, we have developed a new pho- 
tometric redshift method. Comparing the photometric redshifts for 33 galaxies with spectroscopic redshifts 
from VLT/FOS2, we find that the error in the photometric redshifts is less than 0.02{l+z). 

2. Combining multiwave band ground-based and ACS/NIC photometric data, we estimate the photo- 
metric redshifts for galaxies in the UDF. We find the BzK color criterion can be used to select star-forming 
galaxies at 1.4 < z < 2.5 with /iTvcga < 22.0. 

3. Down to i^vega < 22.0, we selected 52 sBzKs in the UDF. Most of the sBzKs have irregular 
morphologies on the ACS images. The logarithm of the number counts of sBzKs increases linearly with the 
K magnitudes, while that of EROs flattens out by K^jcga. ^19. 

4. The two-point angular correlation function of sBzKs and field galaxies were calculated, and the 
clustering amplitudes of sBzKs are found to be about a factor of 6 higher than those of generic galaxies 
at the same magnitude range. A smooth decline in ampUtude with K-hwA magnitude is found, which is 
consistent with results from previous works. 

5. Using the approximate relations from Daddi et al. (2004) and multicolor photometry, we estimated 
the color excess, SFRs and stellar masses of sBzKs. These ifvega < 22 galaxies have a median reddening 
of E{B -V)^ 0.28, an average SFR of ~ 36 Mq yr-\ and a typical stellar mass of ~ 10^°Mq. 

6. The UV criterion for the galaxies at z ^ 2 can select most of sBzKs in the UDF (Xvcga ^ 22.0), but 
it does not work well for the brighter sBzKs sample in the COSMOS (i^vega ~ 19.2). A large fraction of 
massive, highly reddened, actively star-forming galaxies can not be selected by the UV criterion. 
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